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Finite  Element  iaalyiii  of  Slov-ffave 
Schottky  Contact  Printed  Line* 


Abstract 


Extensive  finite  element  analyses  on  MMIC  (monolithic 
microwave  integrated  circuit)  slov~vave  structures  with  both 
localized  and  layered  models  are  presented.  Good  agreement  is 
achieved  for  results  obtained  in  this  paper  with  other 
theoretical  results  and  experiments.  Higher  order  elements  that 
improve  accuracy  are  discussed.  The  comparative  studies  for  the 
Schottky  contact  microstrip  and  the  coplanar  waveguide  with 
localized  and  layered  models  are  presented.  Potential 
applications  of  the  localized  models  to  more  general  and 
practical  slow-wave  circuits  are  also  discussed. 
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Instead,  they  both  have  localized  depletion  regions  on 


semiconductors.  The  effect  of  these  localized  depletion  regions 


has  not  yet  been  discussed.  It  is  plausible  to  contemplate  that 


the  layered  model  commonly  used  for  the  analysis  may  not 


correctly  describe  the  actual  field  distributions  in  the 


structure  with  localized  depletion  regions. 


In  this  paper  the  finite  element  method  based  on  Ez  -Hz 


formulation  [9.11]  is  used  for  a  variety  of  slow-wave 


structures  which  can  be  best  described  by  the  localized  model. 


In  conjunction  with  this  study,  use  of  higher  order  elements 


,e.g.  quadratic  isoparametric  elements, is  discussed  and  results 


are  presented. 


CHAPTER  2  :  the  Localised  Models 


The  localized  depletion  models  for  Schottky  contact 
microstrip  and  CPV  are  illustrated  in  Figure  2.1  and  Figure  2.2 
respectively.  By  setting  the  conductivity  in  Region  III  of  Figure 
2.1  and  that  in  region  IV  of  Figure  2.2  to  zero,  the  localized 
models  will  be  reduced  to  the  conventional  MIS  or  layered  models. 

Obviously,  the  boundaries  of  the  depletion  region  in  actual 
devices  are  not  straight  lines  but  curved.  Such  curved  boundaries 
can  be  found  from  solutions  of  a  static  Poisson's  equation  for  a 
given  bias  condition.  The  present  algorithm  can  then  be  applied 
to  such  geometry.  However,  the  objective  of  the  present  paper  is 
to  study  the  effect  of  localization.  The  essential  feature  can  be 
found  without  solving  the  structures  with  curved  depletion 
boundaries.  Therefore,  the  localized  depletion  region  is  assumed 
to  be  rectangular  in  shape. 

Furthermore,  there  is  also  practical  interest  in  the 
development  of  low  loss  slow-wave  circuit  elements.  Schottky 
contacts  for  both  center  and  ground  strips  of  the  CPW  are  assumed 
under  the  same  DC  bias  conditions.  It  is  found  that  the  loss  in 
this  CPW  structure  can  be  less  than  the  case  where  only  the 
center  conductor  is  a  Schottky  contact  in  the  CPW. 


Layered  model  :  o3»o4-c5*o6"0 
Localized  model  :  a  3-c  5 -c  6-0 


VI  o6,e6 


setiii conductor  layer  o2*E2  II 


1  al,el 

semi-insulating  substrate 


Figure  2.2  Localized  depletion  model  for  Schottky  contact 


coplanar  waveguide. 


CHAPTER  3  :  Derivation  of  the  7EM  Matrix  Equation 

A.  Theory 

The  conventional  Ez-Hz  formulation  that  results  in  a 
homogeneous  coupled  symmetric  matrix  equation  is  adopted  in  this 
paper  [9*11].  In  the  lossless  case*  the  finite  element  method 
with  Ez-Hz  formulation  provides  a  variational  solution.  The 
matrix  equation  can  be  derived  by  a  functional  followed  by  the 
Ritz  approximation  [11]  or  by  weighted  residual  integration 
followed  by  the  Galerkin's  method  [9]  .  It  is  possible  to  obtain 
a  variational  statement  of  the  problem  even  for  a  lossy  structure 
by  dividing  Ez  and  Hz  components  into  real  and  imaginary  parts 
[12].  Alternative ly,  the  3-vector  H-f ormulation  can  be  used 

[13,14].  This  includes  mixed  boundary  conditions  and  cautions 
must  be  exercised  in  treating  the  singularity  at  the  corner  of 
conductor  edge.  The  conventional  Ez-hz  formulation  used  here  does 
not  provide  a  variational  solution  for  the  lossy  system.  The 
method  used  here  forces  the  residual  to  be  zero  by  making  it 
orthogonal  to  each  member  of  a  complete  set  of  the  trial 
functions.  This  method  is  one  manifestation  of  the  method  of 
weighted  residuals,  which  does  not  require  the  existence  of  a 
variational  principle  [15,16],  The  method  of  weighted  residuals 
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uniform  in  z-direction 
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A  uniform  waveguide  discretized 


by  eight-node  quadratic  elements 
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4;  u0  and  £q  are  free  space  permeability  and  permittivity* 
while  ure  and  £j.e  are  relative  permeability  and  relative 
dielectric  constant*  respectively.  Subscript  e  denotes  element 
number  e  associated  with  subdomain  of  the  entire  waveguide 
cross-section  which  contains  Ne^e  subdomains. 

The  equations  (1)  and  (2)  yield  four  equations: 

VtxEz-YU2xEt  -  -ju)ueHt 
VtxEt  -  -juPe^z 
VtxHz-YUzxHt  -  joieeEt 
VtxHt  -  jo)€eEz 

Let  the  inner  product  be  defined  as 

<R,  S>  *  I  R*S  dxdy 

.and  the  test  functions 

P=POz  .  P=0  on  magnetic  wall. 

Q=QUZ  •  Q=0  on  electric  wall. 


(3) 

(4) 

(5) 

(6) 


/iv.v.  «*.  a  .v.  f.v.v.v.  \  -  «■  v 


it  can  be  shown  from  the  equation  (4)  that 


x  Et, 

V 

H 

-j(aue<Hz,  P> 

<7t 

x  Et, 

-► 

p>  - 

f  +  +  +  •> 

|  n  x  Et*Pds  +  <Et, 

Since 

(i) 

(ii) 


(iii) 


n  x  Et  ■  0  on  electric  wall. 

I  n  x  Et*P  =  0,  P=0  on  agnetic  wall 

'n 


(-*•  ■*  +,  !  -*  ■*  ■* 
j  nj  x  Eti*Pds  +  j  n2  x  Et2"Pds 

re  re+l 

-  j  (nj  x  Etj  -  nj  x  Et2),pds  »  0 

re 


• it  ia  concluded  that 


Nele  ^  ^  Nele 

l  <7txEt,  P>  -  l  <EC,  VtxP> 
e=l  e=l 

Similarly,  from  the  equation  (6)  one  obtains 


The  equations  (3)  and  (5)  can  be  reduced  to 


K|  Et  -  -jwue(7tHz)xUz-Y(?cEz) 
Kl  Ht  “  j^e(7tEz)^z-Y(vtHz) 

•where  K|  -  w2ueee+Y2 


(1 

(1 


Combine  the  equations  (8).  (9),  and  (11)  to  obtain 


1  I  A(ure(VtP)*(7cU)+  j^(7cP)x(VtV).Uz]dxdy 

e*l  Re 

*  (1 
Nele  . 

-  I  I  ureU‘P  dxdy 

Similar  to  the  derivation  of  the  equation  (13).  the  equations 
(6),  (7),  and  (10)  can  be  reduced  to 


r  1  n 

1,1  w(^re(7tQ)'(vtv)  -  j^tQJ^t^’Mdxdy 
e»l 

(1 

Nel« 

"  I  |  £rev*Q  dxdy 
e«l JQe 


Note  that  the  normalized  variables  are  introduced. 


The  variable*  U,  V  and  the  teat  functions  P,  Q  are  then  expanded 


in  the  following  way: 


U  (x,y)«  [  U4W1t 
5=1 


V  (x,y)-  l  ViWie 
5=1 


•  where  Vjm  and  Wje  are  shape  functions  at  j-th  node.  N  is  the 
total  number  of  nodes.  0r=0  if  the  r-th  node  is  on  the  magnetic 


wall  and  Ug=0  if  the  s-th  node  is  on  the  electric  wall. 


P  (x.y)  -  [  PiWi 


Q  (x.y)  -  I  QiWie 


.where  WiB(xr.yr)=0  for  i=l,N  if  the  r-th  node  is  on  the  magnetic 


wall,  and  Wie(xg,yg)=0  for  i=l.N  if  the  s-th  node  is  on  the 


electric  wal  1. 


The  final  FEM  homogeneous  matrix  equations  are  obtained  by 


substituting  the  equations  (15)  through  (18)  into  the  equations 


(13)  and  (14). 


Ptk(^U  +  Ptk(2)v  -  0 


9TK(3)y  +  gTK(4)v  -  o 


A'V.vV-.v.’V  ■-,V 


K  NxN  “ 

1  r2_^'7twim) '  (7twjm^  "  yrewimwjm 

,:xay 

] 

V  Ke 

k(2)- 

z  NxN 

L  ^Vtwim  *  9tWje  dxdy 

-NxN 

f  x  VtWjmdxdy  -  [K‘^]T 

K  ^  m 
~  NxN 

f  ["“=!  (vtwie)  *  (7twje)  “  erewiewje 

dxdy 

Oh  e 

The  equation  (20)  can  be  written  in  a  compact  form. 


[PT  qT] 


KfD  K(2) 


k(3)  k(4) 


"u 

V 

J 

. 

-  0 


Since  the  test  function  functions  P  and  Q  are  arbitrary. 
(21)  can  be  reduced  to 


For  a  non-trival  solution  of  X,  det(A)  must  be  zero. 


C.  Application  of  the  Finite  Element  Method 

The  important  steps  in  the  actual  coding  of  the  final  matrix 
equation  obtained  can  he  found  in  [17].  The  hierarchy  of  the 
program  can  accept  any  isoparametric  elements.  In  the  present 
paper*  four-node  bi- linear  and  quadratic  eight-node  elements  are 
used.  The  four-node  element  has  complete  linear  polynomials  along 
its  boundary*  and  the  eight-node  element  has  complete  quadratics. 


15 


Their  differences  in  terms  of  numerical  results  will  br 
discussed. 


By  alternating  Ez  and  Hz  nodal  variables  in  the  column 
vector  X,  the  matrix  A  becomes  a  banded  matrix.  Final  ly,  rows  and 
columns  corresponding  to  the  Dirichlet  boundary  conditions  are 
deleted  such  that  Ez  and  Hz  vanish  at  electric  and  magnetic 
vails,  respectively.  Therefore,  the  matrix  A  has  dimension  MxM. 
where  M  equals  (2N-L).  N  is  the  total  number  of  nodes  and  L 
stands  for  the  sum  of  the  total  number  of  nodes  located  on 
electric  or  magnetic  wills.  The  complex  root  of  the  equation. 
det(A)=0.  is  the  solution  for  the  propagation  constant.  The  real 
and  imaginary  parts  of  the  propagation  constant  correspond  to  the 
attenuation  constant  and  the  slow-wave  factor,  respectively. 


CHAPTER  4  :  Numerical  Results 


A.  Validity  Check 


The  present  FEM  code  was  applied  to  lossless  printed  line 


structures  and  excellent  agreement  obtained  with  the  available 


data  such  as  Fig.2.7  and  Fig. 7. 11  in  [18].  For  the  lossy  layered 


case,  Figure  4.1  and  Figure  4.2  compare  the  results  obtained  by 


the  present  method  as  applied  to  a  layered  structure  and  other 


existing  data  for  MIS  CPV  slow-wave  propagation  both 


theoretically  and  experimentally.  The  small  discrepancy  among 


various  methods  may  be  attributed  to  the  CPW  structure  which  has 


relatively  high  aspect  ratio  (b/a),  i.e.  10,  and  very  thin 


insulating  region.  None  of  the  numerical  solutions  can  find  very 


accurate  answers  over  the  frequency  span  of  interest.  In 


particular,  it  is  well  known  that  the  slender  element,  which  is 


very  long  on  one  side  and  very  short  on  the  other  side,  may  yield 


poor  accuracies  in  the  FEM  analysis  of  a  physical  problem  [19]. 


Owing  to  the  fact  that  either  the  depletion  or  the  lossy 


semiconductor  region  is  extremely  thin,  it  is  inevitable  that 


slender  elements  exist  if  the  use  of  higher  order  and  finer 


elements  is  not  incorporated  in  the  FEM.  Therefore,  quadratic 


8-node  elements  and  total  number  of  96  nodes  are  used  for  this 


„*■*_*»  ^ Vs* » **  * *»**“ *•>*'■  *' - /*  »**  *  *  * 


particular  structure. 


The  case  for  MIS  microstrip  slow-wave  propagation  is  shown 
in  Figures  4.3  and  4.4,  where  the  discretization  based  on  bi¬ 
linear  elements  and  total  number  of  36  nodes  is  sufficient  to 
match  the  data. 

B.  Bi-linear  versus  Quadratic  Elements 

Figures  4.5  and  4.6  compare  computational  results  for  slow- 
wave  propagation  in  MIS  microstrip  model  based  on  both  bi-linear 
and  quadratic  elements.  Both  slow-wave  factor  and  attenuation 
constant  agree  very  well  in  low  substrate  conductivity,  say  less 
than  .1  S/m. 

When  bi-linear  elements  are  used,  the  slow-wave  factor  and 
attenuation  constant  approach  infinity  as  substrate  conductivity 
increases.  Clearly  this  is  non-physical,  because  very  highly 
conducting  substrate  can  be  regarded  as  a  metal  layer  and  the 
slow-wave  factor  should  be  brought  down  when  the  substrate  turns 
into  metal.  In  fact,  in  the  case  of  quadratic  elements,  the 
slow-wave  factor  starts  to  decline  gradually  and  the  attenuation 
constant( loss)  starts  to  increase  again  when  substrate 
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Figure  4.3  Validity  check  of  the  results  obtained 


by  this  paper  and  the  data  in  the 


literature  [9];  the  slow-wave  factor 


versus  frequency. 
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Figure  4.4  Validity  check  of  the  results  obtained 

by  this  paper  and  data  in  the  literature 
[9];  the  attenuation  constant  versus 


frequency. 
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Figure  4.5  Comparisons  of  numerical  results; 

the  slow-wave  factor  versus  substrate 
conductivity  based  on  quadratic  and 
bi-linear  elements. 


(1)  dl-.20a  d2-2ui  2a-.  10m 

(2)  dl-.20M  d2>2ua  2a-.10aa 
l  .  (3)  dl-.25M  d2-l#ja  2i-.1Em 


:  bi-linear 
:  quadratic 
:  quadratic 


10  .1  l 

substrate  conductivity 


Figure  4.6  Comparison  of  numerical  results; 


the  attenuation  constant  versus 


microstrip  substrate  conductivity  based 
on  quadratic  and  bi-linear  elements. 
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conductivity  is  increased  £rom  approximately  500  S/m.  Similar 
comparative  studies  are  performed  for  localized  CPW  models.  The 
same  phenomena  that  both  slow-wave  factor  and  attenuation 
constant  approach  infinity  as  the  substrate  conductivity 
increases  are  observed  when  bi- linear  elements  are  employed. 

Another  type  of  comparative  studies  have  been  performed, 
i.e.  numerical  computations  of  the  slow-wave  factor  and 
attenuation  constant  versus  frequency  by  applying  both  bi- linear 
and  quadratic  elements  on  the  same  CPW  structure.  The  results 
are  illustrated  in  Figures  4.7  and  4.8.  In  these  figures,  curves 
based  on  quadratic  elements  are  identical  to  those  used  in 
Figures  4.1  and  4.2.  The  solutions  obtained  by  means  of  quadratic 
elements  are  apparently  much  closer  to  results  obtained  by  other 
theoretical  methods  and  experiments  than  those  of  linear 
elements.  Additional  studies  on  MIS  microstrip  also  draw  the 
same  conclusion. 

Since  the  accuracy  of  the  quadratic  elements  is  confirmed, 
the  rest  of  the  computations  in  the  paper  are  based  on  quadratic 
elements  except  the  results  shown  in  Figures  4.9  and  4.10. 
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Comparison  of  numerical  results, 
the  slow-wave  factor  versus  frequency 
on  the  same  MIS  CPW  structure  analyzed  in 
Figures  4.1  and  4.2  based  on  quadratic 


and  bi-linear  elements. 
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Figure  4.8  Comparison  of  numerical  results, 

the  attenuation  constant  versus  frequency 
on  the  same  MIS  CPW  structure  analyzed  in 
Figures  4.1  and  4.2  based  on  quadratic 
and  bi-linear  elements. 
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Effect*  of  the  Localized  Depletion  Regions 


on  the  Slov-Vave  Propagation 

It  is  clear  in  Figures  4.9  and  4.10  that  the  MIS  layered 
model  is  a  good  model  for  microstrip  slow-wave  propagation.  This 
conclusion  is  confirmed  again  using  quadratic  elements  with 
improved  accuracy.  The  solutions  obtained  by  layered  and 
localized  models  are  very  close  and  cannot  be  distinguished  by 
plots.  Table  4.1  shows  how  close  the  solutions  are  under 
different  combinations  of  conductivities.  For  the  localized  CPW 
model,  the  results  shown  in  Figures  4.11  and  4.12  are  rather 
interesting.  Two  CPWs  with  various  dimensions  are  syudied  here. 
Structure  (2)  has  about  half  the  gap  width,  i.e.  (b-a),  to 
structure  (1)  and  all  the  rest  of  the  conditions  are  almost  the 
same.  The  results  indicated  taht  the  conducting  Region  IV  has  a 
stronger  influence  on  the  structure  (2).  A  noticeable  transition 
region  exists  where  the  attenuation  constant  (loss)  starts  to 
decline  and  the  slow-wave  factor  changes  abruptly  and  settles 
quickly  as  04  increases.  The  skin  depth,  which  is  inversely 
proportional  to  the  square  root  of  the  conductivity,  is 
relatively  large  in  Region  IV.  approximately  50  mm  for  04  = 
0.1  S/m.  Electromagnetic  field  penetrates  Region  IV  freely  and 
interacts  with  it  in  terms  of  lossy  dielectric  material.  When  34 
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Figure  4.9  Comparison  of  plots  of  the  slow-wave 
factor  versus  frequency  based  on 
localized  and  layered  models  for 
Schottky  contact  microstrips. 
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Figure  4.10  Comparison  of  plots  of  the  attenuation 
constant  versus  frequency  based  on 
localized  and  layered  models  for 
Schottky  contact  microstrips. 
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Figure  4.11  Effects  of  the  localized  depletion 


region  on  CPW's;  the  slow-wave  factor 
versus  the  conductivity  of  the  region 
four(IV)  of  Figure  2.2. 


.01  .1  1  10  100  1000 

conductivity  in  region  IV(=-> . S/m 


Figure  4.12  Effects  of  the  localized  depletion 
region  on  CPW's;  the  attenuation 


constant  versus  the  conductivity  of 
the  region  four(IV)  of  Figure  2.2. 


is  increased  to  10^  S/m  the  skin  depth  is  reduced  to  50  um.  This 
is  very  close  to  the  gap  widths  in  structures  (1)  and  (2).  As  a 
result,  the  highly  conducting  Region  IV  can  be  in  effect  replaced 
by  a  piece  of  thin  and  imperfect  metal  inserted  into  the  gap.  The 
loss  introduced  at  higher  values  of  04  is  much  more  closely 
related  to  skin-effect  ohmic  loss.  Generally  speaking,  the  loss 
mechanism  due  to  Region  IV  in  the  CPW  structure  can  be  separated 
into  dielectric  loss  and  conductor  loss.  The  former  increases  as 
conductivity  increases  and  the  latter  decreases  as  conductivity 
increases.  Common  waveguides  and  transmission  lines  also  exhibit 
similar  loss  behavior  [20]. 

Qualitative  discussions  can  be  readily  extended  to  obtain  an 
equivalent  transmission  line  circuit  representation  of  the 
localized  CPW  in  Figure  4.13.  This  empirical  model  happens  to  be 
a  slight  modification  of  the  analytical  model  of  the  Schottky 
contact  coplanar  line  based  on  semi-empirical  considerations  and 
is  indeed  identical  if  we  combine  Ci  and  C3  into  a  single 
capacitor.  Note  that  in  the  localized  CPW  structure, 
semiconductor  Schottky  contacts  are  under  the  metal  strips.  L0 
and  Rq  represent  the  inductance  and  resistance  due  to  ohmic  loss 
and  skin  loss  per  unit  length,  respectively.  C0  represents  the 
capacitance  per  unit  length  between  center  and  ground  strips  in 


Figure  4.13  Approximate  transmission  line  circuit 
representation  of  the  Schottky  contact 


CPW  with  localized  depletion  regions 
under  metal  strips. 
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Region  VI.  Cj  represents  the  capscitance  per  unit  length  from 
center  strip  to  the  edges  of  Region  IV  and  Region  II.  C3 
represents  the  capacitance  per  unit  length  from  the  ground  strip 
to  the  other  side  of  the  Region  IV  and  the  edge  of  Region  II. 
Parallel  R2  and  C2  represent  the  circuit  contributions  of  Region 
IV  and  part  of  Region  II.  General  discussions  leading  to  the 
determination  of  these  component  values  can  be  found  in  [21.22], 

Figures  4.14  and  4.15  compare  numerical  results  obtained  by 
layered  and  local  ized  model  s  for  the  slov-vave  propagation  on 
CPU's.  The  localized  model  tends  to  show  better  performance  with 
more  linear  slow-wave  characteristics  versus  frequency  and  less 
loss  at  higher  frequency  in  the  structure  under  analysis.  It 
suggests  that  the  structure  with  Schottky  contacts  for  all  metal 
strips.  Figure  2.2.  may  become  a  practical  slow-wave 
configuration. 
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Figure  4.14  Comparison  of  plots  of  the  slow-wave 


factor  versus  frequency  based  on 


localized  and  layered  models  for  CPW’s. 
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Figure  4.15  Comparison  of  plots  of  the  attenuation 
constant  versus  frequency  based  on 
localized  and  layered  models  for  CPW's 


CHAPTER  5  :  Applications 


Along  with  the  study  of  the  effects  of  the  localized 
depletion  models  on  MMIC  slow-wave  circuits  a  general  purpose  FEM 
program  was  developed  for  the  analysis  and  design  of  planar 
waveguides.  Figure  5.1  shows  several  possible  applications  among 
a  broad  class  of  MMIC  waveguides  that  can  be  analyzed  by 
extensions  of  the  present  code.  These  structures  can  not  be 
modeled  accurately  by  an  MIS  layered  model.  The  results  obtained 
from  the  analysis  of  CPW  and  microstrip  MMIC's  with  localized 
models  suggest  that  it  is  possible  to  tune  the  slow-wave  circuit 
to  meet  certain  design  specifications  by  selective  ion 
implantations  on  the  locations  proposed  in  Figure  5.1  with 
appropriate  controls  over  doping  concentration  and  geometry.  A 
certain  mode  of  propagation  may  be  more  subject  to  the  existing 
selective  ion- impl anted  regions  than  the  other  mode  of 
propagation  and  the  additional  ion- implanted  region  can  enhance 
or  reduce  the  loss  associated  with  it  depending  on  what  the 
doping  profile  is. 

The  FEM  code  can  also  be  applied  to  analyze  a  variety  of 
waveguides  with  non-uniform  material  propertities.  These  non¬ 
uniformities  can  be  characterized  by  the  input  data  with  various 
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Figure  5.1  Potential  applications  of  the  extensions  of  the 

localized  depletion  models  to  MMIC  planar  waveguide 
designs,  e.g.  coupled  coplanar  lines  and  coupled 
lines  with  selective  ion  implantations  which  add 
more  degrees  of  freedom  to  design  slow-wave  circuits. 
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material  properties  knovn  at  each  specific  node.  For  example,  an 
optically  controlled  CPW  phase  shifter  proposed  by  Cheung  and 
Neikirk  [23]  has  non-unifom  carrier  concentrations  in  the  semi¬ 
conducting  region  after  optical  injection.  Such  non-uniformities 
can  be  effectively  analyzed  by  the  FEH  field  analysis  program. 


CHAPTER  6  :  Conclusions 


Extensive  computer  simulations  employing  the  FEM  technique 
on  CPW  and  microstrip  MM1C  slow-wave  structures  have  been 
performed  to  study  the  differences  on  slow-wave  propagations 
between  layered  models  commonly  used  in  the  past  and  the 
localized  model  which  is  closer  to  the  physical  situation.  The 
results  show  that  the  layered  model  is  a  good  approximation  for 
the  microstrip  case.  On  the  other  hand,  the  localized  model 
yields  more  accurate  solutions  for  Schottky  CPW  structures.  It 
was  also  pointed  out  that  higher  order  elements  improve  solutions 
for  the  geometry  under  analysis.  The  work  can  be  extended  to 
analyze  and  design  more  complicated  planar  waveguides. 
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layered  models  are  presented.  Potential  applications  of  the  localized 
models  to  more  general  and  practical  slow-wave  circuits  are  also  discussed. 
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